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appearing in this equation is the sum of two terms: a Liouvillian corresponding to a 
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the diffusion constant characterizing the asymptotics of the motion of the particle is 
proven to be strictly positive and hnite. If the disorder in the random potential is so 
large that transport is completely suppressed in the limit where the thermal noise is 
turned off, then the diffusion constant tends to zero proportional to the coupling of 
the particle to the heat bath. 
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I. INTRODUCTION 


In this paper, we study the propagation of a quantum particle through the lattice l/' in 
the presence of a disordered potential landscape and under the influence of thermal noise. 
The dynamics of the particle is described by a Lindblad equation for the time evolution of its 
state, a one-particle density matrix. The Hamiltonian part in the total Lindblad generator 
is given by a random Schrodinger operator (i.e., an Anderson Hamiltonian), the dissipative 
part, which describes the dynamical effects of the thermal noise, by a Lindblad operator that 
couples the motion of the particle to the degrees of freedom — afterwards “traced out” — of a 
heat bath in thermal equilibrium at a positive temperature If the coupling of the particle 
to the heat bath is turned off, its dynamics is generated by a standard random Schrodinger 
operator. It is well known that strong disorder in the potential landscape may completely 
suppress coherent transport of the quantum particle — the particle gets stuck near a well of 
the random potential. This is the phenomenon of Anderson localization. It is natural to ask 
whether localization survives when the particle is coupled to the heat bath. The main result 
established in this paper (see Theorem 11.11 below) says that an arbitrarily small amount of 
thermal noise suffices to destroy Anderson localization. The particle then exhibits a diffusive 
motion (“quantum Brownian motion”). We will further show that, under the assumption of 
complete localization in the absence of thermal noise, the diffusion constant tends to zero, 
as the coupling to the heat bath is turned off. In contrast, in the absence of disorder, the 
diffusion constant diverges, as the coupling between the particle and the heat bath is turned 
off, because, in this limit, transport of the particle is governed by ballistic motion. 

This paper is an adaptation to the present context of work of the second author on 
dynamical randomness^ and makes use of similar mathematical techniques. 

To be concrete, we consider a quantum particle hopping on the lattice whose pure states 
are given by wave functions, 'ipt G ^^(Z*^), evolving according to the Schrodinger equation 

dt li’t) = \ipt) 

where t denotes time. In this equation, is an Anderson (random Schrodinger) Hamiltonian 
of the form 

= u I^)( 2 /I + ^5Z uj{x) |x) (a:| (1.1) 

\x-y\ = l X 
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where u is the hopping amplitude, a; is a random potential, with {uj{x)}^(zzd independent iden¬ 
tically distributed (i.i.d.) random variables, and A is a constant characterizing the strength 
of the disorder. It is often assumed that the distribution of the variables u{x) has a bounded 
density. In this paper, we only need to assume that if the disorder is large all states are 
localized in the sense that the spectrum of is dense pure-point, and the correspond¬ 
ing eigenfunctions are localized. This is known if the distribution of the variables u{x) is 
bounded. Assuming that the ratio “/a is sufficiently small, there then exists a constant /i > 0 
such that if 'ipt=o has finite support, e.g., {x\'ipt=o) = then 

E(|(x|^/.t)|') < Ce-^l"l (1.2) 

where E (■) denotes averaging with respect to the random potential uj, and C* is a finite 
constant (depending on the choice of 'ipt=o)', see Frohlich and Spencer - and Aizenman and 
Molchanov-. See also Ref. 0 for a recent survey with an estimate of the critical value of “/a. 

In concrete situations of physics, a quantum particle will usually interact, at least weakly, 
with degrees of freedom describing a surrounding medium, e.g., with the phonons corre¬ 
sponding to quantized vibrations of a crystal lattice. If these degrees of freedom are excited 
it is conceivable that, even at large disorder, the particle exhibits dissipative transport — 
diffusive motion, or motion with friction — due to its interactions with the medium. This is 
the phenomenon studied in this paper. 

A physical theory of dissipative electron transport in conducting materials was devel¬ 
oped a long time ago; see, e.g., Mott-. More recently the mathematical formalism for the 
quantum Markov approximation in aperiodic and disordered media was studied by Spehner 
and Bellissard-, although quantum diffusion was not established in that work. A simple 
one-particle model exhibiting quantum diffusion was studied in Ref. 0, without a quantum 
Markov approximation. A formalism for the study of dissipative transport in disordered 
semi-conductors, within a quantum Markov approximation but incorporating effects of the 
Pauli principle on a gas of non-interacting electrons at positive density, was developed by 
Androulakis et al.-. The purpose of the study undertaken in Ref. 3 was to gain some under¬ 
standing of the role of thermal noise in the quantum Hall effect. In comparison, the goals of 
the present paper are more modest. Assuming that the density of particles (moving, e.g., in 
a conduction band of a semi-conductor) is so small that the one-particle approximation can 
be justified, we propose to analyze the interplay between randomness and thermal noise in 
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the transport of a single qnantnm particle. We plan to analyze systems at positive density 
in fnture work. 

Next, we introdnce some notation and describe the model studied in this paper. The state 
of an open system consisting of a single quantum particle that interacts with a dynamical 
environment (medium) is typically a mixed state, which is described by a one-particle density 
matrix, pt (with t denoting time). This one-particle description is obtained by tracing out 
all degrees of freedom of the environment. In the kinetic (Markovian) limit, the evolution of 
Pt is described by a a Lindblad equation of the form 

dipt = -i[H^,pt]+ gCpt, (1.3) 

where is a Hamiltonian, in the following the one introduced in eq. o, g is a constant 
that measures the strength of the thermal noise, and the operator £ is a Lindblad generator, 
which we specify below. 

In our context, a one-particle density matrix, p, is a positive, trace-class operator on 
In the a:-space representation, the operator p has matrix elements {x\p\y), with 
{x, y) G X It is convenient to introduce the variables 

X = x + y, ^ = x-y. (1.4) 

Then X G Z'^ and ^ G Z'^, with 

X ± ^ G 2Z'^ . 

From now on, we will write p{X,^) for the matrix element (x| p \y) of p, where it is under¬ 
stood that (X,^) and {x,y) are related by eq. fll.4j) . Consistent with eq. fll.4j) we adopt the 
convention that p(X, ,^) = 0 if X or X — do not belong to 2Z'^. 

To dehne the action of the Lindblad operator C on the density matrix p we introduce two 
operators, G (for “gain”) and L (for “loss”), dehned as follows: 

(Gp)(A',?) = (1.6) 

and 

(Lp)(.Y,f) = ^r(0,i,-f)p(.Y,),) (1.6) 

where the “gain kernel” r : Z'^ x Z*^ —)■ C satishes certain properties specihed below. More 
generally, the operators G and L might also act on the variable X G Z'^ by convolution 
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(which is imposed by the requirement of translation invariance), but this possibility will not 
be considered in the present paper. 

The Lindblad operator C is expressed in terms of G and L by 


(£p)(X,0 = {Gp){X,i)-{Lp){X,i) 


(1.7) 


so that the Lindblad equation fll.Sp for pt, incorporating hopping, disorder and dissipation, 
takes the form 


dtpt{X, Cj = - in ^ [pt{X + e, ^ + e) - pt{X + e, ^ - e)] 

|e|=l 


— iA 


u - \ — u 


p{x,i) 


( 1 , 8 ) 


+ 9 ')) - '■(0,1) - 0] 9(V, I)) . 


The gain kernel r is usually assumed to have certain physically natural properties, such 
as “detailed balance” for the energy. These properties are discussed in Al below. Here we 
only note that, for our analysis to proceed, r must satisfy the following 

Assumption I.l. The gain kernel r : x —)■ C satisfies: 

h) = 0 unless ^ + p E 22f; 


2. r the Fourier transform of a non-negative measure p, 

r{f,v) = [ e‘«'P"‘'^'‘^/i(dp,dq), 

where denotes the d-torus = [0, 27r)'^; 

3. r{f^, 0) — r(0, —0 = 0 for each G Z"^; and 
4- for any f G £^(Z'^) we have that 

Re ^ ^ [r(0,r7-O-r-(00] 0(00(0 > c^|0(Or 

5. r is the kernel of a bounded operator on 0(Z'^), i.e., 


E 






for some G < oo and any f G £^(Z'^); 


(1.9) 


(I.IO) 
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Remark. 1) Item 1 is required for the gain operator {Gp){X,^) to satisfy the constraint of 
vanishing whenever X ± ^ 27/“. Items 1 through 5 are discussed in §I Al below. 

2) Since r is the kernel of a bounded operator on it follows that £ is a bounded map (from 
—)■ on each fiber over X G Without loss, we scale r and g such that C has 

norm < 1, i.e., 

For simplicity, we further assume that the Lindblad generator C respects certain lattice 
symmetries: 


Assumption 1.2. 

i = 1,..., d, 


h) invariant under inversion of coordinates. That is, for each 


r{Rii,Riv) = f (?.>)) 


(M2) 


where 


RiiXli ■ ■ ■ 1 ■ ■ ■ I ^d) 


(Xi, 


-Xi, 




is invariant under permutation of coordinates. That is, for any permutation a 
of{l,...,d}, 

r{T„f,TaV) = r{^,v) (1.13) 


where 

£7(^1? • • • : Xffj (.^(7(1)7 • • • : 


From Assumption 11.21 and eq. fll.Sp it follows that the processes (t, X, pt{T„X,T„^), 
(t,X,^) I—)■ pt{RiX,Ri^) and (t,X,^) 1 —)■ pt{X,^) all have the same distribution. 

The following theorems are the main results proved in this paper. 


Theorem I.l. Let pt he a solution of eq. Ci with initial condition pt=o = |0) (0|, and 
suppose that r satisfies Assumptions \I.l\ and U.Si If the particle interacts with the heat bath, 
i.e., g > 0, then the diffusion constant 


Dij := XiXjE{{x\ pt\x)) = ^ X^Xj-E (pt(X, 0)) (1.14) 


cgZ'^ 


xe2Z'i 


exists and satisfies 


Dij = D5ij 


(1.15) 
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d 

with 0 < D = < oo. 

i=l 

Let us now consider the behavior of the diffusion constant in the limit where the coupling 
with the heat bath tends to zero, g ^ 0. We are interested in the behavior of the diffusion 
constant D as a function of g, so we will henceforth write D{g). In particular, we would like 
to estimate the size of D{g) if the hopping u and the disorder strength are such that the 
strong localization result described in eq. flI.2D holds. In fact, we will not use the full strength 
of eq. flI.2D . Instead, what is required below is simply a uniform bound on the second moment 
of the position, viz. 

I |2 

:= sup ^-^E ^|(a;| |0)|^j < cxD . (1-16) 


Theorem 1.2. Suppose eq. fll.lbp holds. 


Then 


D{g) = ^g + o{g) , (1.17) 

with 0 < A < (1 + Vc)f'^. (Here c is the constant appearing in item 5 of assum,ption I/, il ) 

By way of contrast, in the absence of disorder (A = 0 in eq. fll.lD ) we have the following 
elementary result. 

Theorem 1.3. //A = 0 then 

Dig) = C- (1.18) 

g 

with 0 < C < Yc. 

We observe that the diffusion constant diverges, as (7 —)■ 0; the reason being that, without 
disorder and without dissipation, the motion of the particle is ballistic. 


A. Properties of the gain kernel 

Let dq denote normalized Haar measure on the d-torus T"*, dq = 1. The measure fi 
appearing in eq. flL9j) need not be absolutely continuous with respect to the product measure 
dp X dq on x Nonetheless, we will write (using distributional notation) 

h(dp,dq) = f(p,q)dpdq 
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where r is a (possibly singular) positive distribution on x — see eq. fll.29p below. In 


terms of r, the gain and loss operators have the following expressions 


(Gp'^)(X,p) := [ 

JT<i 


p'*'(V,q)r(p, q)dq 


(1.19) 


and 


(Lp^)(X,p) := f f(q,p)dq p^{X,p) 


( 1 . 20 ) 


Here is the Wigner transform of p, which is the Fourier transform of p{X, in the f 
variable: 


p"'(A',p) - 5^e‘i> V(.Y,?) 


( 1 . 21 ) 




The variable p is the quasi-momentum of the particle; it ranges over the d-torus T'^, which 
is the Brillouin zone corresponding to 7/“. Although it need not be positive, the function 


p^(X, p) is interpreted as a density for the distribution of the particle position and quasi¬ 
momentum in phase space x 


When acting on the Wigner transform the Lindblad operator has the form 


(£p^)(X,p) := (Gp^)(X,p)-(Lp^)(X,p) 


/ [^(p, q) - ^(q, p)] dq ( 1 . 22 ) 

jfd 


Thus, in the hber over W, C acts as the generator of a hopping process in the quasi¬ 
momentum: r(p, q)dq is the rate at which the quasi-momentum of the particle jumps from 
q to p. To be able to interpret C as the generator of a stochastic process, we need r to be a 
positive measure. Thus we take r to be of positive type — item 2 of Assumption 11.11 
Item 3 of Assumption II.II is equivalent to 



(1.23) 


Thus the generator of the jump process may just as well be written as 



(1.24) 


for / G It follows that £1 = 0 where 1 is the identity function l(p) = 1 on the torus, 

i.e., Haar measure dp is an invariant measure for the momentum jump process. In fact, by 







item 4 of Assumption lI.il Haar measure is the unique invariant probability measure. Indeed, 
eq. fll.lOl) is readily seen to be equivalent to the following estimate: 


-2Re(/, £/)^2(Td) = [ ^(P,q)dpdq|/(p)-/(q)l^ 

>c[ dpdq|/(p) -/(q)l^ (1.25) 

for every / G C{T‘^ x T'^). Eq. fll.25p is a “spectral gap” condition for the generator of the 
jump process in quasi-momentum space. Assuming that this condition holds amounts to 
requiring exponentially fast mixing of the jump process. 

Item 3 of Assumption II. II implies that the operators 


GfiO = and L /(0 = '^ r{0,r] - ^)f{r]) 

are bounded on A natural condition on the jump process which gives the bound is 

to require the total rate of jumping out of any hxed quasi-momentum p to be uniformly 
bounded, i.e.. 


M : = 


:= sup / r(q, p)dq < oo . 

P 


( 1 , 26 ) 


Lemma I.l. If the measure r(p,q)dpdq satisfies eqs. (II.23P and (11.2611 then G and L are 
hounded operators on 


Remark. That is, eqs. (11.23^ and (11.26^ are sufficient for item 3 of Assumption 11.11 


Proof. Taking Fourier transforms we hnd that 


G' 0 (p) = / f(p,q) 0 (q)dq and L 0 (p) = / f(q,p)dq 0 (p) 
where 0(p) = is the Fourier transform of 0. By Plancherel’s Theorem, 

IIL^IIj < MUW., , 

with M as in eq. (II.26p . To bound G0, we note that 


G0(p) 


< M 


Ifd 


0 (q) f(p,q)dq 


by the Cauchy-Schwarz inequality, since by eq. (II.23P 


M = sup / r(p, q)dq. 

peT'^ 
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Integrating over p and applying Plancherel’s Theorem again, we find that 


The results we have arrived at, so far, are summarized in the following 

Proposition 1.2. Let r(p, q)dpdq be a positive measure on x T'^ that is periodic with 
period tt under joint translations of p and q, i.e., 


f(p + 7rn, q + 7rn) = f(p, q) 


(1.27) 


If r is defined by eq. (PI) and eqs. fll.23|l . fll.25p and flL26|) hold for r, then r satisfies As¬ 
sumption [TJ 


Eq. (11.271) is required for item 1 of Assumption lI.il namely that r{^,'r]) = 0 if ^ + 77 ^ 2Z^, 
to hold. In this regard, it is useful to note that 

p"’(X,p+n^) = 

for n G since in the sum over ^ we have X — ^ E 2lfi. Thus, is either periodic or anti- 
periodic with period vr in each component of the quasi-momentum, depending on whether 
the corresponding component of X is even or odd. Eq. (II.27^ guarantees that Gp^ and Lp^ 
also have this property. 

Lindblad operators are often assumed to satisfy a “detailed balance condition” at the 
temperature, /3~^, of the heat bath. For a free particle detailed balance can be formulated 
as follows. Let e(p) denote the energy of a freely moving particle with quasi-momentum p 
— i.e., 

/ d 

£(p) = 2u \ d — cos(pj) 

\ i=i 

for the isotropic nearest neighbor hopping considered here. Detailed balance at temperature 
is the condition that 

^(P) q) = p) . (1.28) 

For example, for a Lindblad generator describing the interaction of the particle with a thermal 
reservoir of non-interacting bosons with an even dispersion law cj(k) = a;(—k) at temperature 
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(3 one has that 


^(p,q) = ^(p-q) 




^-^ij(p-q) 


+ 


1 — e“/^^(P“q) 


^(5(£(p) - e(q) + u(p - q)) 


(1.29) 


where S{-) denotes a Dirac delta fnnction and the “form factor” F is a non-negative, even 
fnnction. 

In the present context detailed balance, as in eq. fll.281) . is only consistent with item 4 of 
Assnmption[LT]if (3 = 0 (infinite temperatnre). At /3 = 0, detailed balance is just symmetry of 
the gain kernel, r(p, q) = r(q, p). However, symmetry is stronger than eq. flI.23H and detailed 
balance and symmetry, as such, play no role in our analysis. An explicit, and symmetric, 
example of a gain kernel satisfying our requirements is given by a suitable limit of eq. flI.29D 
as /3 —)■ 0, i.e.. 


(p,q) = ^(p-q) 


5(£(p)-£( q)-n;(p-q)) + 5(e(p) - £(q) + n;(p - q)) 


(1.30) 


For a quantum particle whose dynamics (in the absence of thermal noise) is governed by a 
random Schrodinger operator of the kind studied in this paper, detailed balance is a somewhat 
awkward condition. If the disorder is large it is natural to neglect the kinetic energy term 
of the particle in fll.28p . i.e., to set £ = 0. In this case, our Assumption II. II is acceptable. A 
(not very natural) example of a gain kernel to which our analysis would apply, with (3 > 0, 
can be found in Ref. 1^. But these matters ought to be studied more thoroughly. 


II. THE DIFFUSION CONSTANT FOR LINDBLAD DYNAMICS WITH 
DISORDER 


In this section we prove Theorem 11.11 To begin with, we note that by Assumption 11.21 

XiXjE{pt{X,0)) = 0, 


for all t. Indeed, let 


Pi(A,0 = pMX,R,0 
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where Ri denotes inversion of the coordinate. Because pt and pt have the same distribution, 
one has that 

X X 

However, the dehnition of 'pt implies that 

5^X,X,E(pi(X,0)) = -5^X,X,E(pi(X,0)) 

X X 

if z 7 ^ j. Likewise, by permutation symmetry, all diagonal matrix elements are the same 

5^XfE(p,(V0)) = 5 ^VJE(p,(V0)) 

X X 

z, j = 1,... ,d. Thus the off-diagonal elements of the diffusion matrix Dij dehned by eq. 
(11.141) vanish, and the diagonal elements are all equal, provided they are well dehned. Thus, 
to analyze the diagonal elements of the diffusion matrix, it suffices to consider Di i. 

Note that pt is a random variable depending on the disorder conhguration a; G fl. In what 
follows, we will emphasize this by writing pt{X,^,u). The initial condition is the density 
matrix po given by 

Po(X,e,n;) = 5o(^)5o(Ol- • 

We introduce the Hilbert space 

•H = {T G X X ff) I T(X,^,a;) = 0 if X±^ ^2Z'^} . 

Finite-group-velocity estimates for the propagation of a quantum particle on the lattice 
show that 

5^e"^l^l|A(X,e,n;)|'<e^-^ (H.l) 

Here m and Cm are positive, hnite constants; see Lemma lA.ll In particular pt G "H, for each 
t > 0. 

Translation of the system by a lattice vector a is given by a unitary map on R dehned by 

= T(X-2a,e,r_,a;) 

where TaUj{x) = u{x — a). We dehne the following generalized Fourier transform on R: 

XT(e,n;,k) := + 2a, C, r.o;) . 
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This Fourier transform leads to a direct-integral decomposition of "H, with hbers isomorphic 
to l-i := X n). Given T G Tf, we let Tk(^,a;) := ca, k). We then have that 

Tk G 'H, for almost every k, and 



Lemma II.1. Suppose that pt{X,^,u) solves eq. fll.Sp . Then 

dtPk;t = + gCkj Pk;t (II-2) 

where A]^ = uTk 4- XV, with 

|e| = l 

= (w(0-^(0))/(^,^^), and 
5. £k/(^,a;) = ^ [r{^,v) - r{0,p - 0] f{p,Tn^u). 

Proof. We define operators T, V and £ on "H by 


Tp{X,^,u;) := p{X + e,f + e,uj) - p{X + e,f - e,uj) 

|e| = l 


Vp{X,^,uj) := (^u -UJ p{X,^,uj) 

and 

Cp{X,i,uj) := '^[r{i,p)-r{t],r]-f)]p{X,p,uj) 
so that the equation of motion for pt reads 


dtPt = -iuTpt - iXVpt -4 gCpt ■ 


After Fourier transformation, this becomes 

dtAa = -'^'^Tp^^^-iXVp^.^ +gCp^.^ . 
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Straightforward computations yield 


[^(^ + 2a + e,^ + e,raw)-^'(^ + 2a + e,^-e,raw)] 


|e| = l 


E 

|e|=l 


y^e + 2a + e,^ + e,raa;) - - e + 2(a + e),^ - e,Taw) 


[$k(^ + e,u)- - e, T-eUj) 


|e|=l 


E + e, o;) - e + e, Teu) 

|e|=l 


and 


l/d'k(^,w) = E® '^'''{raUj{^ + a)-TaUj{a))^{^ + 2a,^,TaUj) 

a 

= (^^(0-^(0))5^e"'^'“d^(^ + 2a,C,raa;) = (a;(0 - a;(0))§k(^, w) . 

a 

Furthermore, 

£d'k(^,w) = y^e"'*"'“y^[r(^,r7)-r(0,77-O]^(^ + 2a,r/,raCu) 

a 77 

7 ] a 


V 


□ 


We set 


Qk ■ i-^fc 9 ^k • 

Lemma [II. II has the following corollary. 


(11.3) 


Lemma II.2. 


jY^XIMpAXM) = -S|.( 3 o®l, e-«‘ 3 „®l>, 3 |^_„ 


(11.4) 


X 


Remark. Here 


5o®l((r,a;) = 


1 if X = 0 , 

0 otherwise. 
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Proof. First note that, by eq. flll.l|) . we are justified in interchanging differentiation and 
summation to write 


.YfE (p,(A', 0, (a(A', 0, i^)) 


X 


X 


k=0 


By shift invariance of the distribution of u, it follows that 


5^X2EU(X,0,a;)) = - e^'^'^E (p,(2a;, 0, r.o;)) 


X 




k=0 


Since pq{X,^,cj) = (5o(-^)^o(0) have that Pk,o(^,i^) = ^o(0- Lemma [11. II then implies 
that pk,t = e“*®''5o 0 1. This completes the proof. □ 


To compute the right hand side of eq. (111.411 . we begin by noting that 

-dl (5o0l, 

= -2 [ ds rdr(ho0l, 

Jo Jo 

+ [ ds{So®l, dlg^\^^^e-^^°6o(^l) 


(11.5) 


= -2 / ds I dr(dk,gl 5o 0 1, e 9fci^klk=o <^0 0 l) 


'0 JO 
rt 


k=0 


+ / ds(5o0l, ^fci^klk=o^o'^ 1) 


where we have used that 


0 1 = e“*^ohn 0 1 = 0 1 . 


( 11 . 6 ) 


Eq. (III.6P follows from the following facts 


1 . ^ 0^0 0 1 — 0 , 

2. £k<^o 0 1 = 0, and 

3. 4(5o0 1 = 0, 


which may be verified by explicit computation. (In fact, (1) and (3) are quite general — 
they follow from the conservation of quantum probabilities. However, (2) requires Item 1 of 
Assumption 11.111 
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Since ® 1 = ® 1 = 0, for all k, we find that 

® 1 = iudkiTkSo (g) 1 = M [e‘^^5-ei 0 1 - 0 l] 

5fci^k'^o0l = -m<9fciTk5o0l = -M [e'*'^h-ei 0 1 - 0 l] 

and 

^fci^k^o 0 1 = m [e‘^^5_ei 0 1 + e~'^^de^ 0 l] 
where ei is the hrst basis vector. Using these computations in eq. flll.Sj) . we conclude that 

((5o0l, e“*^‘‘<5o 0 = 2u^ / ds /" dr (0, 

where 

0 = (hei - 5-ei) 0 1 • 

Using Lemma III.21 it now follows from the Tauberian theorems formulated in Ref. U, 
Chapter XIII, that 


'0 ^'O 


lim ^ Vx2E(pRX,0,a;)) = lim 2 ^ 2 ! f ds f dr U, e'^^V) 

t^oo At ■++ t^oo t Jo Jo ^ ^ 

POO nt 

= lim2M^77 / dfe“*^ / dr (0, e“^®° 


X 


r ?—>0 


'0 


'0 

= lim 2 m^ [ dre~^^((j), = \im2u^ ((j), -—( 

provided the limit on the right side exists. 

This limit can be shown to exist using a straightforward “Feshbach argument.” Recall 
that ^ 0^0 0 1 = Go^o 0 1 = 0. It follows that 

Y = {fe X I (<^0 0 1,/) = 0} 

is invariant under Go- Since 0 G we need only consider the restriction of Go to 
to compute the resolvent matrix element of the resolvent we wish to control. We further 
decompose the subspace 'R“‘“ as follows: 


where 


Y = 


e 




f{x, Cj) = 0 if X 7 ^ 0. i = {5o 0 / I E (/(cj)) = 0} 
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and 


ni = [f{x,u)en f{o,u) = oj . (11.7) 

Let Qq denote the orthogonal projection of onto "Hq and set Qi := 1 — Qo- We have the 
following block-matrix form for the action of Qq on 


fo\n^ - 


0 iuQoToQi 

iuQiToQo QiGoQi 


Since (j) G Tff, the Schnr Complement formnla yields 


’ h + ^o 


V + QiGoQi + ^QiTqQqToQi 


Now 


Re QiGoQi = gQi^GQi > cgQi 
where c is the constant appearing in eq. fll.iop . Thns 


1 

Tj + QiGoQi + ^QiTqQoTqQi 



1 

< — < oo 

eg 


uniformly in rj. 

It follows that (0, (?7 -f Go)~^G) remains bonnded as r/ —)■ 0; and it remains to show that 
it has a non-zero limit. To this end, let 


n := orthogonal projection of T-L^ onto kernel of QiTqQqTqQi 
= orthogonal projection of onto kernel of QqTqQi ■ 


Then we have (see Lemma IB.1( below) 

lim-;-- = If———LI, 

^^0 V + QiGoQi + f QiToQoToQi RGoR 

where the limit is in the weak operator topology and the operator fl^on is bonndedly- 

invertible on the range of fl, becanse 

Refl^on = flTTI > cgll . 


By an explicit compntation, we see that given f{x,oj) G 


QofoQif 


(0,a;) 


- fie,reUj)) ■ 

\e\=l 
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It follows that 


QoToQl (<5ei - 5-ei; 

Hence 0 = {Se^ — 5-ei) ® 1 G ranH, and 


1 = 0 . 


lim 

r ]^0 


v + Go 




Thus, 


D := hm - ((a;| Pi lx)) = 


HQoU 


Note that D < and D > 0, since it is a limit of positive quantities. Moreover, 


D = Re D = 2gu^Re ( d), C- 


UGoR 


> 2gu^c 


RQoR 


> 


Agu^c 


since ||0||^ = 2. This completes the proof of Theorem 11.11 


III. DIFFUSION WITHOUT DISORDER — THEOREM 0 

In the last section, we have not made use of the disorder in our estimates. Indeed, the 
result also holds when A = 0. In fact, for A = 0, the calculations become much simpler, 
because we have that Mo = Tq = 0, and hence 


D = 2u (cj), 

\ gC 

By eq. fll.lOp . (0, ) is bounded above by ‘^/c. Theorem 11.31 follows. 


IV. PERTURBATION THEORY FOR D IN THE LARGE DISORDER 
REGIME — THEOREM n 


Localization in the form of Eq. fll.lbp implies that 


D(0) = lim -^Ixil^E 

t—^OO t ^ ^ 


{x\ e 


| 0 )| 


= 0 . 


In order to compare the diffusion constant at p = 0 with the one at p > 0, it is convenient 
to introduce a cut-off version of D{g), namely 

poo 

D{g,g) := df e"^* V jxil'E ((x| p* |x)) (IV.l) 

^0 
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where pt satisfies eq. fll.3j) . The quantity D{g,ri) can be thought of as the diffusion constant 
on a time scale of Y??. By Lemma [II.21 


D{g,v) := dl (So^l, 


T] 


Af + gC + g 


-(5n(8) 1 


/c=0 


Following the calculations of the previous section, and using that 

u(j) = 5A:ii-4fc|fc=o(^o <8 1 and <8 1 = 0, 

we see that 


D{g,g) = 2u^((j), 


iA + g^ + v 


(IV.2) 


and 


D{g) = \im D{g;g) . 

rj^O 


For g = 0, 


D{0,g) = 

Jo 


^IxipE (|(x|e-‘^“Y0)'^ 


dt, 


hence 


D{0,g) < gi'^ 

with i is as in eq. (1I.16F As D{0,g) is real and Aq is Hermitian, eq. (1IV.2P implies that 


D{Q,g) = ReD{0,g) = 2gu^ 


i^o + g 


Thus 


and so 


i^o + V 


g 

2m2 




^|x|2E(|(x|e-^^“Y0)'^ 


dt 


i^o + g 


2 n2 

< - 

“ 2m2 


Because Aq is a Hermitian operator, 


1 . 

' - f 

1 

iA + g'*' 

Jm 

it + g 


where is the spectral measure of Aq associated to <p. By monotone convergence, 

2 n2 


= lim 


iAo + g 


< 


2u^d 
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Thus, -0 = exists in H, and 


< 


2m2 


Note that 


'ijj = lim 


r?—>0 —iwAo + Tj 

where ■0 is the complex conjugate of ^fJ. (This is true, because 0 is a real function and is 
a real symmetric operator.) 

Translated back to the Schrodinger operator picture, we have proven the existence of the 
limit 


:= \imri 

rj—>-0 




(|(a;|e-‘*^“ |0)| 


and obviously io < 1. Note that 
Returning to (/ > 0, we have that 


= D{0,ri) -2u^g {(j), ^ - C, ^ 


12 £2 

I 2u^d' 


dt 


i^o + V + V 


+ 2nV(0, T 




-C- 


i^o + V i^o + 5''^ + h ^*^0 + V 

Since C = QiCQi, where Qi is the orthogonal projection onto the subspace introduced 
in (1II.7I1 . we have that 

C __ C = C __£ . 

ivAo + gC + g iQiA-oQi + gC + g + jjQiToQoToQi 

Taking the limit g ^ 0, as in the proof of Theorem lI.il we hnd that 

1 


D{g) = -2u g ( ip, 


C - gCTl- 


n£ 


tp 


inA-oif + 

where If is the projection onto the kernel of QqTqQi, as above. In particular, we obtain the 
upper bound 

\D{g)\ < 2u^g (l + -'] W^jf = (l + -\ilg 


where we have used that ||£|| < 1 (by convention). Furthermore, we have the lower bound 

Au^gc Au^gc 


Dig) > 


11^0 


Uo + gCr 
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obtained at the end of the proof of Theorem II.ll We conclnde that 


D{g) = 0{g), as 0 . 


(IV.3) 


In fact, we can take the limit 


A := l.mkM 


3 to g 

Let IIo denote the projection onto the kernel of Il^on. If IIo = 0 then 


(IV.4) 


9 




> 0 


ill^on + gliCIi 

in the weak operator topology; see Lemma fB.ll It seems likely that Hq = 0, bnt we are not 
aware of a proof. Fnrthermore we do not need to resolve this issne, since if IIq 7 ^ 0 then 

g 




■n„ 


ill^QlI + ^n/TI IIq/TIq 

in the weak operator topology; see Lemma IB.ll In any case, the limit in eq. flIV.411 exists, 
and 


A = —( Ip, 


£-£nn 


1 


no>^iLio 


-nn£ 


ip 


where we take no(no/TIo) ^IIq = 0 if IIq = 0. Note that we have the estimates 


4c 


IIA"^ - 


< A < 1 + - L 
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Appendix A: Finite gronp velocity estimates for the Lindblad eqnation 


Lemma A.l. Let po be a given density matrix such that 

>1 = supe’“l^l <«). 

If, for t > 0, the density matrices pt satisfy eg. (O. with initial condition po, then 

supe^W 

with Cm = 4:de^u + g. 

Proof. Let Bm denote the Banach space of fnnctions F : x —)■ C such that 

Ill’ll™ - supfee^^inA'.Ol'j < 

Thus Po G Bm by assumption. Let = iX (cj — ta (^^)) and let 

Note that |lTi(X, ^)| = |pt(X, ^)|. In particular, Wq E Bm- 

The evolution of Wt is governed by the non-autonomous equation 

dtWt = GtWt , 

with time-dependent generator 

Gt = + gC)e-'^^^ , 

where the “kinetic hopping operator” T is 

Tp(X,0 = 5]|p(A + e,f + e)-p(X + e,^-e)| . 

|e|=l 

Because cf and C hber over X, one easily computes that 


||ei*'^£e-i*'^F| 


supe^l^l I 


Y, (>•(?. >)) - r(0,,, - 0) ,,) 


< ll^ll, 
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by the normalization eq. fll.llj) . On the other hand, the hopping term e^^'^Te is bounded 
by 


||e'*'^Te“'*'^F| 


= supe 

X 


m\X\ ( ^ 


+ e, ^ + e) - + e,^-e) 

|e|=l 


< 2supe’"W5^ O^|F(A' + e,0r 

|e|=l V t 


X 


< 2 


sup 


X 


|e|=l 


II^L < 4de™||F||, 


Thus 


\\GtF\\^ < Me-u + g 


and the result follows. 


□ 


Appendix B: A limiting principle for resolvents 

Lemma B.l. Let LL be a Hilbert space. Let A be a normal operator on H and B a bounded 
operator on H, with Re A > 0 and Re R > c > 0. 

1. If kei A = {0}, then 

lim (0, (XA + By^fj) = 0 

for any 0, -0 G "H. 

2. If ker A ^ {0}, then 


lim (0, (aa + r)"VL = (n0, (nRn)^^n0) 

A^oo ' ft \ ' 


ran 11 


for any 0,0 G "H, where If = projection onto the kernel of A. 


Proof. Let 0 G be given and let hx = (AA + B) ^0. Note that ||hA|| < c ^ ||0||, since 


c||hA||^ < Re(hA, 0) < II^aI 


23 









This and the identity \Ahx = -tp — Bhx imply 


IAIMAaII < (l + c-'IIBII) WII . 

Thus (J — n)hA converges weakly to zero. If ker A = {0}, so 11 = 0, then this completes the 
proof. On the other hand, if 11 7 ^ 0, then it commutes with A, since A is normal. Thus 

UBhx = . 

Since (/ — n)hA converges weakly to 0 and Ili?!! is boundedly invertible on ran 11, it follows 
that n/iA converges weakly to □ 
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